Phase transition in hierarchy model of Bonabeau et al 
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The model of Bonabeau explains the emergence of social hierarchies from 
the memory of fights in an initially egalitarian society. Introducing a feed- 
back from the social inequality into the probability to win a fight, we find 
a sharp transition between egalitarian society at low population density and 
hierarchical society at high population density. 
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How do social hierarchies emerge in (human) society ? Bonabeau et al 
[1] explained it in the tradition of statistical physics, i.e. as random. They 
claimed a phase transition between egalitarian society at low population den- 
sities and hierarchical society at high population densities. Sousa and Stauf- 
fer [2] found this transition to be an artifact of the non-stationarity caused 
by their assumptions; and when the rule was changed to give stationarity, 
the sharp transition vanished. Now another simple rule change is shown to 
recover a phase transition, though at a different population density. 

In the model on Bonabeau et al [1] people initially are all equal and are 
distributed randomly on a square lattice such that no two people occupy 
the same place, filling a fraction p of the lattice sites. They then diffuse 
randomly, and whenever one person i wants to move onto the site of person 
k a fight between these two people occurs. If % wins, then % and k exchange 
their sites on the lattice; if instead k wins, none of the two moves. (The 
population density is thus kept constant during the whole evolution, and is 
used as an extrinsic parameter determining the final behaviour of the society: 
egalitarian or hierarchical.) The probability q for i to win is given by 

q = l/(l + exp(r][h(k) - h(i)]) (1) 

where rj is a free parameter and h(j) counts the weighted number of victories, 
minus the weighted number of losses, of person j. This weight is initially 
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Order parameter versus density; L = 1 00 (lines), 300 (+), 1 000 (x) and 5000 (squares); t=1 0000 
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Figure 1: Variation of "order parameter" a versus population density p, for 
10 4 iterations in L x L lattices, L = 100, 300, 1000 and 5000. 

unity and then decreases by ten percent at each iteration (= one attempt 
to move on average all people through random sequential updating). Thus 
the memory of past fights fades away according to an exponential rate, i.e. a 
finite-time memory. For low population density the memory is lost between 
the rare collisions; for high density the memories of the now frequent fights 
add up. Simulations [2] gave a smooth increase of the amount a of inequality 
with increasing density p. This amount of inequality is measured by 

a = (< q 2 > - < q > 2 fl 2 (2) 

and initially is zero if all q are 1/2 (since all h — 0). 

To get a sharp transition, instead of a smooth variation, in <j(p), we now 
replace 77 in Eq.(l) by a, i.e. we introduce a feedback of the order parameter 
a to the rules of the fights [3]: 

q = 1/(1 + exp(a[h(k)-h(i)}) . (3) 
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In reality this means that losing or winning has more influence on the future 
in an hierarchical than in an egalitarian society. (For the first ten iterations, 
the factor a is omitted since then no inequalities had yet the chance to build 
up.) Fig.l shows our results for lattices of sizes 100 x 100. . .5000 x 5000, 
with a rather clear transition at a critical point p c ~ 0.32. Fig.2 shows better 
data (averages over typically 64 lattices) close to the transition; it seems the 
order parameter jumps from zero to about 1/4 when p increases above p c . 

We thank Suzana and Paulo Murilo de Oliveira for discussion, G. Weis- 
buch for help during a computer breakdown, and the Jiilich supercomputer 
center for time on their Cray-T3E. 



Averages over (mostly) 64 lattices; L = 100 (lines), 300 (+) and 1000 (x) 
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Figure 2: As Fig.l, but restricted to p ~ p c and with better statistics. Note 
the systematic downward shift of p c {L) with increasing lattice size L. 
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